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Abstract 

We construct integrable generalised models in a systematic way exploring different repre- 
sentations of the gl{N) algebra. The models are then interpreted in the context of atomic and 
molecular physics, most of them related to different types of Bose-Einstein condensates. The 
spectrum of the models is given through the analytical Bethe ansatz method. 
We further extend these results to the case of the superalgebra gl{M\N), providing in this way 
models which also include fermions. 
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1 Introduction 



Exactly solvable models are a fascinating issue that continue to attract considerable interest in physics 
and mathematics. Altough the integrability of quantum systems is usually restricted to one dimen- 
sion, there are many reasons that turn this study relevant for physical applications. It serves as a test 
for computer analysis and analytical methods for realistic systems to which, until now, only numer- 
ical calculations and perturbative methods may be applied. In addition, a nontrivial solvable model 
reveals the essence of the phenomena under consideration. For instance, many concepts established 
in critical phenomena were inspired by the exact solution of the Ising model. From the experimental 
point of view, there are some real spin-1 compounds (e.g. NENC, NDPK, or NBYC etc. . . ) [1] and 
strong coupling ladder compounds (such as (5IAP)2CuBr4-2H20, or Cu2(C5Hi2N2)2Cl4, etc..) [2] 
that can be perfectly well described by integrable models. The necessity of using exactly solvable 
models has been also demonstrated through experimental research on aluminium grains at nanoscale 
level [3]. 

A significant aspect of integrable systems is its interdisciplinary character, i.e., they can be found 
in different areas of physics. The Ising and the Heisenberg models [4] in statistical mechanics, the 
t-J and Hubbard models [5] in condensed matter physics, the nonlinear cr-model [6] in quantum 
field theory, the interacting Boson model [7] in nuclear physics and more recently the two-site Bose- 
Hubbard model [8] in atomic and molecular physics are just some representative examples of the high 
impact and potentiality of these systems. Let us remark also the emergence of integrable systems 
in high energy physics, more particularly in gauge theories [9-11] (for a recent review, see [12]), or 
string theory, through the recent analysis in super- Yang-Mill theories, see e.g. [13,14]. 

Therefore, new exactly solvable models are highly welcome and constitute the main focus of the 
present article. In particular we will concentrate on the construction of integrable generalised models 
in atomic and molecular physics, most of them related to Bose-Einstein condensates (BECs). 

The phenomenon of Bose-Einstein condensation, while predicted long ago [15], is currently one of 
the most active fields in physics, responsible for many of new perspectives on the potential applica- 
tions of quantum systems. Since the early experimental realizations of BECs using ultracold dilute 
alkali gases [16], intense efforts have been devoted to the study of new properties of BEC. In recent 
years the creation of a molecular BEC form an atomic BEC has been obtained by different tech- 
niques [17]. The field was further broadened by the achievement of quantum degeneracy in ultracold 
fermionic gases [18]. These achievements could lead to new scientific investigations that includes 
coherent atomic lasers, quantum chemistry, the quantum gas with anisotropic dipolar interactions, 
quantum information, atomtronics, among many others. 

In this context, it is natural to expect that exactly solvable models in the BEC scenario may 
be of relevance, providing some physical insights [19]. Our main purpose here is to employ the 
integrable systems machinery in its full power, i.e., exploring all possible types of representations of 
some algebra (we consider, in particular the gl{N) algebra) to enlarge the family of known exactly 
solvable models in atomic and molecular physics with the aim that potentially new relevant models 
emerge. Using this machinery some existing models in the BEC scenario, such as the two-site Bose- 
Hubbard model [8] will be restored as well as new ones will be obtained. Remarkably, a two-coupled 
BEC model with a field, a two-coupled BEC-model with different types of atoms and a three-coupled 
BEC model, among others, will be introduced in this general framework. It is worth to mention here 
that the popular " BEC-transistor" in atomotronics uses a BEC in a triple well [20]. The models are 
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then solved by means of Bethe ansatz methods. 

Our paper is organized as follows: in section 2 we briefly review the general setting of integrable 
systems and fix notation. In section 3 we present our general approach and also discuss the different 
representations of the gl{N) algebra that will be adopted. Section 4 is devoted to the discussion of the 
different physical models we can get using this construction. In section 5 the Bethe ansatz equations 
of the models are derived. We extend these results to the case of the superalgebra gl{M\N) in section 
6 and some applications of this formalism, i.e. models which also include fermions are presented in 
section 7. Section 8 is devoted to some concluding remarks. 



2 Generalities 

2.1 Monodromy and transfer matrices 

We remind here the general setting used in the context of integrable spin chains and more generally 
in QISM [21-24]. One starts with a so-called 'algebraic' monodromy matrix T{u), which is a x 
matrix taking values in an algebra A: 

N 

T{u) = Tij{u) Eij with Tij{u) G A ; Eij e End{C^) 

where Eij are the N x N elementary matrices (with 1 at position and elsewhere). In the main 
part of the present paper, the algebra A will be the Yangian of gl{N), Y{N), but we will also study 
sup er- Yangians . 

The monodromy matrix obeys the so-called FRT relation [25]: 

Ruiu - v) T,{u) T2iv) = T^iv) T,{u) R^iu - v) (2.1) 

where we have used the standard auxiliary space notation, e.g. Ti[u) = T{u) ® Iat, where Itv is the 
identity matrix and Ri2{u) is the R matrix of A. The i?-matrix obeys the Yang-Baxer equation 

Ri2{ui - U2) Risiui - Us) i?23(w2 - ^^3) = -^23(^^2 - Us) Risiui - %) -^12(^^1 - U2) ■ (2.2) 
It will also be unitary 

Rl2{Ui - U2) R2l{u2 - Ui) = /(^i, M2) I ® I (2.3) 

where f{ui,U2) is some known function. 

It is well-known that one can produce monodromy matrix for several sites by applying the co- 
product AT{u) = T{u) ® T{u) = T^^^ (u) T^'^^ (u) where the superscript labels in which copy of the 
algebra T{u) acts. More generally, one can consider 

T{u)=T^^\u) T^^\u) . . . T^^\u) 

as a monodromy matrix. This is the base of spin chain models, the copies of the algebra defining 
(upon representation) the L quantum spaces (sites) of the chain, see e.g. [26-29] and references 
therein. In what follows, we will call T["](ii) the elementary monodromy matrices, the product of all 
these elementary matrices providing the 'real' monodromy matrix. 
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Let us stress that to fix a physical model, one has to represent the monodromy matrix, i.e. assign 
to each of the T["l(tt) a representation of the algebra A: changing the representations will lead to 
different physical models. 

Once we have a (represented) monodromy matrix T{u), then fl2.ll) ensures that the transfer matrix 
t{u) = trT{u) obeys 

[tiu),tiv)] = 0. (2.4) 

Expanding t{u) in the variable u leads to commuting integral of motions, one of them being the 
chosen Hamiltonian of the system. For instance, in spin chain models, one takes this Hamiltonian 
to be 7i = t(0)~^ ^'(0), while e.g. in the Bose-Hubbard model, it is simply t(0). The other quantities 
constructed from the transfer matrix just produce conserved quantities. If the number of such 
(independent) conserved quantities is sufficiently large, the system is said to be integrable. 

2.2 Automorphisms of the monodromy matrices 

We present some automorphisms of the relation (12. ip that will be of some use in the physical models 
we will study. 

The first automorphism is built on the transposition: starting from a monodromy matrix T{u), 
it is easy to show that 

T{u) T\-u) i.e. Tij{u) Tji{-u) 

where the transposition is done in the auxiliary space. The proof relies on the unitary relation (12.31) 
of the i?-matrix. We will call this automorphism the sign-transposition. 
Another automorphism is the conjugation by a constant matrix 

T{u) MT{u)M-^ with M e End{C^) 

which is a consequence of the invariance of the i?-matrix 

Ml Ms Ri2{x) = Ruix) Ml Ms 

We will call this automorphism a conjugation. 

A particular case of conjugation is the dilatation automorphism 

T{u) —* aT{u) with a G C 

2.3 Hermiticity 

The elementary monodromy matrices we will consider below are always hermitian 

{t{u))^ = T{u) i.e. Tl{u)=T,,{u) (2.5) 

This implies that 

(T{u)f = (rW(M)TPl(w)y = T[2l(M)TW(ix) (2.6) 
so that the total monodromy matrix is not hermitian. However, using cyclicity of the trace, we get: 

t\u) = fr(^T[2](n)TW(M)) = tr {t^^\u) T^''\u)^ = t{u) (2.7) 
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Thus, the transfer matrix is hermitian. This property is vahd only when L = 2 (and L = 1), cychcity 
being not sufficient to get hermiticity of the transfer matrix as soon as L > 3. For this reason, we 
will focus below on the case L = 2, hence ensuring hermitian Hamiltonians. We shall see that even 
with this restriction, we will get most of the models used in the BEC context, as well as new ones. 

3 Bosonic gl{N) models 

We present here the general approach we use, focusing on the case of the Yangian Y{N) = Y{gl{N)) 
[30]. Other cases are presented in the sections below. The i?-matrix we consider takes the form 
[4,31,32]: 

i?i2(a;) = I®I--Pi2, (3.1) 

X 

where 

N 

is the permutation operator. 

3.1 Conserved quantities 

Since L = 2, the monodromy and transfer matrices have an expansion 

2 

Tkliu) = XI ^"^ir^ Tl^^=uJkki and e C (3.2) 

n=0 

t{u) = t2U^ + tiu + to with t2 e C (3.3) 

For BEC models, one generally uses to = ^(0) as an Hamiltonian, while ti and ^2 correspond to 
integrals of motion. It is easy to see that t2 is just a number, but from the explicit form of the 
i?-matrix, one can get other conserved quantities. 

Indeed starting from the relation (12.11) and projecting on the basis elements Eij ^ E^i in the 
auxiliary spaces, one gets: 

[T.^iu) , ni{v)] = (T,^iu)Tuiv) - nj{v)Taiu)^ . (3.4) 

Then, taking i = j, summing on j, and looking at the coefficient of v, one gets 

[t{u) , T,W] = {uj, - u,) (4°) + u T,«) (3.5) 

This proves in particular that the quantities 

h = Ti}^ , Wk = l,...,N (3.6) 

commute with the transfer matrix and are in involution. Thus, they generate integrals of motions. 

Let us remark that, following the value of the Uk parameters, one could get more conserved 
quantities (through the tI:^\ k ^ I, generators), but they will not form an abelian subalgebra: they 
will generate a symmetry algebra for the model. 
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3.2 Spin chain monodromy matrices 

The elementary monodromy matrices for the Yangian can be speciahzed to gl{N) monodromy matri- 
ces using the so-called evaluation map. This amounts to take these elementary monodromy matrices 
to be of the form 



N 



or in matricial form 



L{u) = Lij{u) Eij with L 



/ u + eii ei2 ei3 

621 U + 622 623 



u6ii + e 



(3.7) 



Liu) 



eiN 
e2N 



\ 



\ em 



eN2 



eN-l,N 
eN,N-l U + CNN J 



(3.1 



Here, eij are gl{N) unrepresented generators obeying 



[cjj , Cfc/] — 6jk Cii — 6ii ekj 



(3.9) 



It is easy to show that L{u) obey the relation fl2.ip with the i?-matrix (13.11) . Hermiticity of L{u) is 
ensured by e|j = e^j. Moreover, one can use the Yangian shift automorphism u —>■ u + w to get extra 
free parameters. Thus, 

T{u) = (u + wi) L[2] (u + W2) 
leads to hermitian integrable models with transfer matrix 

N N 



t{u) = ^(n + wi + efj) ® {u + W2 + efj) + ^ efl 



J2] 



(3.10) 



In the context of spin chain models, the parameters Wj are called inhomogeneity parameters. 

As already stated, it is the choice of a gl{N) representation for each of the sites that will determine 
the physical model one wishes to work on. When the representations are highest weight finite 
dimensional ones, it leads to spin chains models. They have been extensively studied and we just 
repeat here well-known facts to illustrate the techniques we shall use with different representations. 

For instance, one can take the fundamental representation of gl{N) 



TT e 



I] 1 



i,j = 1,...,N 



for both elementary monodromy matrices, leading to a (well-known and somehow trivial) two-site 
spin chain. Specifying furthermore to the case of gl{2), one recovers the Pauli matrices 



0-4 



7r(^ei2j 

leading to an Hamiltonian 



TT 



(621) 



0"_ 



TTleii - 622) 



<7z 



7r(eii + 622) 



h 



H = t{0) = (7+ + 



1 

az + - 
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Of course, one could choose another representation, for instance, for gl(2), take the spin 1 represen- 
tation 










7r(ei2)=| 1 1=5+ ; 7r(e2i) = | 100|=5_ ; = \ | (3.11) 


7r(en)= I i 1 =1(5, + I3) ; ^{622)= ( i 1 =hl,-S,) (3.12) 
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leading to the Hamiltonian 



H = t(0) = S+®S^ + S^®S+ + \s,(^S, + \ 

Zi Zt 

Note that in both cases, the parameters Wn do not play any role because the spin chain is too simple. 
For the same reason, the above Hamiltonians coincide with t(0)~^ ^'(0). 

3.3 Oscillator monodromy matrices 

The above framework can be generalized to other (infinite dimensional) representations of gl{N). 



3.3.1 Bosonic and fermionic representations of gl{N) 

We start with couples of creation/annihilation operators (ai,at), with commutation relations 

[tti, a]] = //j 5ij ; [a\, a]] = [a^, aj] = . (3.13) 
From these relations, it is straightforward to check that C{u) defined by 

N 

C{u) = ^ Cijiu) Eij with Cij{u) = ^iU 6ij + — a] aj (3-14) 

obey the relations 

[Cij{u) , Cm{v)] = ^ ^ ^ {Ckj{u)Cii{v) - Ckj{v)Cii{u)j . (3.15) 

It is equivalent to 

Ri2{u - v) Ci{u) C2{v) = C2{v) Ci{u) Ri2{u - v) with Ri2{x) = I®I--Pi2 (3.16) 

X 

so that 

T{u) = C^'\u + wi)C^^\u + W2) 

provides an integrable model. 

In fact, this calculation is valid for an arbitrary number of site^, and it just corresponds to a 
choice of (infinite dimensional) gl{N) representation 

7r(eij) = alttj , ij = 1, ...,N 
*0f course, we will potentially lose hermiticity of the transfer matrix when L > 3. 
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for the elementary monodromy matrices i\3.7\i . The highest weight is the Fock space vacuum |0 >, but 
the representation is reducible and is an infinite sum of finite dimensional representations with fixed 
'particle number' A/" = ^i- ^^^^ call the corresponding monodromy matrix an 'homogeneous 

oscillator monodromy matrix'. 

Focusing on hermitian elementary matrices, one is led to 



and 



The last equation imposes 

qj = qoe'^' with qo,9je 
The parameter go is irrelevant for C{u), and since 



e*^-* aj and e 



is an invariance of the algebra, one can restrict to the case 

Cij{u) = jJiU Sij + a| ttj with /ij G 



or in matricial form 



C(u) 



t 



jj2U + n2 al as 



\ 



(3.17) 



with rij 



(3.18) 



In general, one takes the values fij = 1 to get canonical commutation relations. Then, C{u) has a 
leading term (in u) which is just the identity matrix, in accordance with the definition of the Yangian 
matrix L{u). 

Let us remark that there exists also a fermionic gl{N) representation where now aj and aj are 



fermionic operators: 



C{Au) = ^iU Sij + al aj with fij G 



(3.19) 



In that case, the representation is finite dimensional, since the oscillators now obey the supplementary 



relations (a^)^ = 



This possibility will be used to produce some fermionic models. 



3.3.2 Inhomogeneous oscillator monodromy matrices 

The above calculation is valid whatever the values of the numbers /ij are. In particular, one can take 
the value fij = 0, for some j & J C [1, A^]: the corresponding £.{u) matrix will still obey (12.1 p with the 
i?-matrix (13.11) . This particular value jij = allows a (scalar) representation aj = G C and aj = a*, 
j G J, of the oscillator algebra. For obvious reason, we will call these operators 'constant oscillators', 
and 'inhomogeneous oscillator monodromy matrix' the corresponding elementary monodromy matrix. 
To distinguish it from the homogeneous one, when needed, we will denote it as A(m) instead of C{u). 

Hence, A(m) is an x matrix built on p = A^ — | J| couples {aj, aj), p being independent from 
A^ (provided it is smaller than A^). Then, T{u) will lead to a transfer matrix based pi +p2 oscillators, 
where pn is the number of oscillators in T^"\u) = A^'^\u), n = 1,2. 
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Let us remark that since we have taken the hmit nj — for j G J, the elementary monodromy 
matrices do not start with Itv, but rather with a non-invertible diagonal matrix. In that sense, we 
are not in the Yangian context anymore. However, since the relation (12. ip is still obeyed with the 
i?-matrix (13. ip . this does not affect the relation (12. 4p . so we are still in the framework of hermitian 
integrable models. The underlying algebraic structure, which is very close to the Yangian, was 
studied in [33] and is called 'truncated Yangians'. Keeping in mind this restriction, we will lose keep 
writing that we are in the Yangian context. 

Note also that the limit fij (and aj,a] constant) can be taken at the very end of the 
calculations. Hence, we can consider a general C{u) matrix, keeping in mind that, to get a A(ti) 
elementary monodromy matrix, and depending on the model one wishes to study, some of the aj, 
operators will be in fact complex numbers aj, a*, and the corresponding fij = 0. For other oscillators, 
one chooses in general = ±1. 

Finally, let us remark that for the fermionic gl{N) representation, since now the aj and aj 
operators obey the supplementary relations (a^)^ = = (a])^, it is not possible to take them as 
non-vanishing constants. Taking all the constants to be zero leads to trivial models, hence, fermions 
are excluded from A{u) when dealing with Y{N). Fortunately, we will see below that one can recover 
them when studying models based on super- Yangians (see section [6]). 

3.3.3 Automorphisms of oscillator algebra 

An automorphism will be used to produce new terms in the Hamiltonians. It exists only for the 
bosonic algebra, and consists in a shift by a constant: 



We will call this automorphism a shift of the oscillator algebra. It can be used to produce boundary 
terms in the different models. 

3.4 gl{N) transfer matrices 

We have seen that we have essentially two types of elementary matrices at our disposal, the matrices 
L{u) and C{u), so that one gets three types of transfer matrix: 



{a,a\fi) 



{a + a,a 



+ a*,/i) , a G C for bosons 



(3.20) 



t{u) = tr {L^^\u + wi) L^'^^u + W2)) 



N N N 



~ n E (4' + 4?) + E 4' 4?' + E (-^ 4? + m e!?) 



(3.21) 




t{u) = tr 



N N N 



~ u^{fii Elf +ni) + ^ Elfa^ai + ^ (^2 Elf + wi rn) 



(3.22) 




t{u) = tr {C^^\u + wi) C^^\u + W2)) 



N N N 




) 



(3.23) 
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where the ~ sign means equahty modulo polynomials in u with constant coefficients. As a notation, 
we have introduced = 7r„(eij), n = 1,2, the representation of the gl{N) generators in L^^^{u) and 

called aj, a], /ij (resp. bj, bj, uj) the oscillator algebras in £[-'^1(m) (resp. in £.^^^{u)); Uaj (resp. nf,j) are 
the corresponding number operators. 

The ffist transfer matrix is just a two-site spin chain, where the first sites carries the 'spin' tti 
of gl{N), and the second site the 'spin' These models (and their generalization to an arbitrary 
number of sites) have been studied for a long time, and we will not consider them here. The two 
other transfer matrices lead to several physical models, depending on the choices of: 

• The gl{N) algebra one considers (i.e. the choice of A^) 

• The gl{N) representation in the L{u) part 

• The characteristic (bosonic or fermionic) of the oscillators in the C{u) part 

• The number of 'constant oscillators' in the A{u) part(s) 

• The values (specially zero or not) of the parameters corresponding to these 'constant oscillators'. 

• The use or not of the different automorphisms. 

The next section is devoted to the presentation of the different physical models one can get from 
these choices. 

4 Examples of gl{N) BEC models 

For simplicity, we now focus on the case where the parameters involved in the elementary matrices 
are real, and normalize the commutator of the oscillators to 1. More general Hamiltonians (still 
hermitian) can be obtained keeping complex parameters, as detailed above. 

4.1 Models based on two by two matrices 

The possible elementary monodromy matrices take the form 

L{u) = 2 1 C[u) = t 1 4.1 

= \ (5a ^^^)=( /a u-n) ^'''^ 

plus possibly the use of shift automorphisms. We have used the sign-transposition and dilation 
automorphisms 

A{u) -A\-u) 

plus a redefinition of the parameter to make A(m) similar to A{u). 

Apart from the spin chain model presented in section 13. 2[ one gets 5 different models. The 
conserved quantities Ij have the general form: 

h = fJ'bi rial + fJ-ai ribi and I2 = Hb2 + /^a2 'nb2 (4.3) 
where (ai, 02) refer to the ffist elementary monodromy matrix, and (61, 62) to the second one. 
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4.1.1 A spin-boson model 



We consider t{u) = trL{u + Wi) A{u + ^2)- Up to irrelevant constant terms, it takes the form 

t{u) = Sz + n) + S;, + n) + (a S+a + a S^a^) + n 

leading to Hamiltonian H = t{0) with conserved quantity 

I = n+^S, (4.4) 

Above (a) denotes the single- mode field creation (annihilation) operator, Sz,S± the atomic in- 
version, rising and lowering operators and Wi is the transition frequency. This model describes 
the interaction of a two-level atom with a single- mode radiation field. It was derived in [34], [35] 
using different methods and it reduces to the Jaynes-Cummings Hamiltonian on resonance and in 
the rotating- wave approximation in the absence of the last term [36]. Despite of its simplicity, this 
model has been a source of insight into a better comprehension of the nuances of the interaction 
between light and matter. It is important to remark that a Jaynes-Cummings model interaction can 
be experimentally realized in cavity-QED setups and also, as an effective interaction in laser cooled 
trapped ions [37]. 

Again, the shift automorphism produces boundary term 

Hbound = f3(^^S,{a^ + a) + a{S+ + S-) + wi{a^ + a)^ (4.5) 

to the transfer matrix. The conserved quantity is modified to 

/' = n+^S, + l3ia^ + a) (4.6) 



4.1.2 Generalised spin-boson model 

We consider t{u) = trL{u + Wi) C{u + W2)- Up to irrelevant constant terms, it takes the form 

t{u) = u{ni + 77-2) + '^Szirii - 712) + S+a\a2 + S_a\ai + Wi{ni + 722) 
leading to Hamiltonian 

H = t{0) = ^S:,{ni - 722) + S+a\a2 + S'_a|ai + Wi{ni + 712) (4.7) 
with conserved quantities 

h=ni + ^S, ; I2 = n2- ^S, (4.8) 

Above the oscillators aj, j = 1,2 denote two radiation fields (two photons, for example) interacting 
with a two-level atom. The atom-field interacting term could be interpreted as a scattering of two 
fields with a two-level atom. Here we mention that if linearly polarised light is used, it is possible to 
have the same transition frequency wi (see, for example [38]). 
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When the oscillators are bosonic, one can use the shift automorphism to add a boundary term: 
Hbound = ai (^^Sz{a{ + ai + ai) + S+a2 + S^al + wi{a\ + ai) j 

+ a2(^^Sz{al + a2 + 02) + S+ai + S-a\ + wi{al + a2)^ + aia2{S+ + S-) (4.9) 
to the transfer matrix. The conserved quantities then become 

I[ = ni + ^S, + ai{a{ + ai) (4.10) 
I2 = ^2 - + a2 (4 + ^2) (4.11) 

4.1.3 Simple heteroatomic-molecular BEC model 

We consider t{u) = trC^^^{u + Wi) C^'^\u + W2), with aj, aj, j = 1, 2 for £W and bj, bj, j = 1,2 for £[^1 
Up to irrelevant constant terms, the transfer matrix takes the form 

t{u) ~ u{ni + ^2) + Uai Ubi + na2 nb2 + a|4^i«2 + 4616201 +Winb + W2 Ua 
where we have introduced the notation 

rial = a\ai ; ribi = blbi ; na2 = a\a2 ; nb2 = 4^2 (4.12) 
ni = Uai + Ubi ; n2 = na2 + nb2 (4-13) 
na = rial + na2 ; rib = ribi + nb2 ■ (4.14) 

It leads to Hamiltonian H = t(0) with conserved quantities 

/i = ni ; l2 = n2 (4.15) 

This Hamiltonian is a particular case of the one presented in the section 14.1.4^ and we postpone the 
physical discussion to this section. 

Finally, using the quantum determinant (see appendix), one can also show from C2 given in 
appendix that naUb is also conserved, so that one finally gets as conserved quantities 

ni = rial + ribi ; n2 = na2 + rib2 ; ria = Uai + na2 ; rib = ribi + nb2 
three of them being independent. 

4.1.4 Heteroatomic-molecular BEC model 

We consider t{u) = trC{u + wi) A{u + W2), with a, and b, b^ for C{u) and c, for A(n). 
Up to irrelevant constant terms, it takes the form 

t{u) = u{na + Tic) + nanc + 0^ nb + (3 {b^c^a + a'bc) +winc + W2 ria 

It leads to Hamiltonian t(0) with conserved quantities 

h = ria + rib ] h = ria + (4.16) 
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It is trivial to check that this Hamihonian (adding terms If and /|) corresponds to the heteroatomic- 
molecular Bose-Einstein condensate model [39] 

Ti. = Uaa nl + Ubb nl + Ucc nl + Uab "-a rib + Uac ria ric + Ubc rib ric 

+ ^ana + ^ib rib + picric + ^ [a) he + 1)^0) , (4-17) 

for the particular choice of the couplings \Uaa = Ubb = U^c = Uq + Ui, Uab = 4[/o, Ubc = 0, 

Uac = 4:Ui + l, Ha = W2, ^b = /?^ /^c = Wi, Q = (3. 

In this context, the parameters Uij describe S-wave scattering, //j are external potentials and 
is the amplitude for interconversion of atoms and molecules. One gets a three-mode Hamiltonian 
describing a Bose-Einstein condensate with two distinct species of atoms, denoted h and c, which can 
combine to produce a molecule a [39]. The total atom number (Ji + I2) and the imbalance between 
the atomic modes {I1 — I2) are conserved quantities. A detailed classical and quantum analysis of this 
model reveals unexpected scenarios, such as the emergence of quantum phases when the imbalance 
is zero [40]. 

4.1.5 The two-site Bose-Hubbard model 

We consider t{u) = trA^^^{u + wi) A[^1(m -|- W2), with a, for A'^^ and 6, for A'^l 
Up to irrelevant constant terms, it takes the form 

t{u) = u{na + rib) +rianb + uj {b^a + a)b) +winb + W2 Ua 

with w = aj3. It leads to Hamiltonian t(0) with conserved quantity 

I = na + nb. (4.18) 

It is easy to verify that by combining the conserved quantity / with the Hamiltonian t(0) and choosing 
properly the coupling constants, we arrive atQ 

H = ^{na-nbf-^{na-nb)-^{a^b + b^a). (4.19) 
o Z Z 

This is the two-site Bose-Hubbard model, also known as the canonical Josephson Hamiltonian [8]. It 
describes the tunneling between two single particle states or modes (a and 6), which can be separated 
spatially (two wells) or internally (two different internal quantum numbers). The parameter K 
corresponds to the atom-atom interaction, A/i is the external potential and £j is the coupling for the 
tunneling. Despite of its apparent simplicity, this model predicts the existence of a threshold coupling 
between a delocalised and self-trapped phase [41], [42], in qualitative agreement with experiments [43]. 

4.2 Models based on three by three matrices 

The number of possibilities increases very fast, we present here only some cases that we found 
physically relevant. The interested reader can easily compute the other models using the techniques 
we have described. 

^More specifically, we have H = c{I^ — 4t(0)), where the following identification has been done: K/8 = c ; 
(A^)/2 = Acw2 = -4cwi ; £j/2 = Acw. 
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4.2.1 Two-coupled BEC model with a single-mode field 

We consider the elementary monodromy matrices 

(u + Ua h a^a^ \ / u + rif. Pi /?2 

6t and A[21(m) = [3ic /3j Pi P2 
a3a asb al / \ P2C P1P2 Pi 

T{u) = A^^\u + wi) A^^\u + W2) (4.20) 

Then, the transfer matrix reads 

t{u) ~ (ric + na)u + Pi rib + 'nanc + Wi + "-a + "3 P2 (a^ c + a) 

+Pi (aU c + a) + «3 Pi P2 {b^ + b) (4.21) 

with conserved quantity 

I = na + Uc- 

This leads to the HamiltoniaiJil 

H = HaUa + ficric + w'^ Hh + (tt^ c + o) + w {tt^ b c + c' b^ o) + w (b^ + b) (4.22) 

which describes two wells {A and C) with atoms interacting with a single-mode field. Above and 
b denote the single field creation and annihilation operators, respectively and and a (c^ and c) the 
creation and annihilation operators for a particle in the well A (C). The total number of particles 
n = Ha + Tic, where Ua = a^a and ric = c^c is conserved. The parameter Q is the coupling for 
the tunneling between the two wells, w is the radiation field frequency and /ic are the external 
potentials. Here there are basically two mechanisms that allow the atoms to trap between the wells: 
i) the tunneling effect, which is related with the depth of the wells; ii) a tunneling which occurs 
mediated by a single-mode field b. 

When w = 0, one has a supplementary relation 

[t{u),nb] = 0, 

and this Hamiltonian reduces to the two-site Bose-Hubbard fl4.1.5p . When f2 = we get 

[t{u),nb -nc\ = 

and this Hamiltonian reduces to (14.1.41) 

4.2.2 Three coupled BEC model 

We consider the elementary monodromy matrices 

/ u + ni a\a2 Pza\ \ / -Pi P1P2 Pia\ 

AW(w) = a\ai u + n2 Pza\\ and k'^'^^u) = \ p^p^ -pl p^al 
\ p3 ai p3 a2 Pl I \ Pl as /?2 03 u- 

T{u) = M^\u + wi) M^\u + W2) (4.23) 



^Thc following identification has been done: Pi ~ w ; fia ~ W2 ; /ic = lui ; = azP2- 
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where in A[^1(m) we have used the sign-transposition and dilatation automorphisms, and a redefinition 
of the parameters Pi and P2. Then, the transfer matrix reads 

t{u) ~ Pi P2 {a\a2 + a|ai) + P3 Pi (ajas + alai) + Ps P2 {ala^ + 0302) - Pini - Pi n2 - Pi 

It leads to the Hamiltonian t{0) with conserved quantity / = 711+77,2 + ^13. This corresponds to 
a special three coupled BEC model with asymmetric tunneling and external potentials. To our 
knowledge, this is the first integrable quantum model describing a three well systeni§. We remark 
here that there has been recently an increasing interest in the study of three-well systems (trimers) 
for a variety of reasons, such as 

i) its possible application in the construction of a BEC-transistor [20]; 

a) it is the simplest model which provides a bridge between the double- well and the multi-well 

systems [44], [45]; 

in) recent achievements in the experimental field, in particular the control promised by microtraps 
[46] suggest the realization of the trimer to be at hand [47] . 

It is easy to check that by combining the conserved quantity / with the Hamiltonian t{0) and 
choosing properly the coupling constants, we arrive at the following Hamiltoniarfl 

H = 1^2 (a2^^i + ala2 + + 0302) + (a|a3 + OgQi) + /i77i + /i773 + yU2 ^^2, (4.24) 

which describes an array of three coupled wells, which will be referred to as the left (1), middle (2) 
and right (3) wells, respectively. Above, Q (resp. Q2) denote the tunneling of atoms between the left 
and the right wells (resp. the left-middle tunnelling and the middle-right tunneling), while fi2 and 
H are the external potentials. Obviously, adapting the choice of parameters, one can also treat the 
case where the left-middle and middle-right tunnellings are different. 

For the particular case where Q approaches to zero this model reduces to the asymmetric open 
trimer model in the absence of the interatomic scattering and a large external potential fi2 [47] . We 
observe here that if we consider this model in its full generality [47], with Q2 and fi2 as adjustable 
parameters, a BEC-transistor [20] can be derived. 

In the symmetric limit ^2 = ^2 = /i2 in fl4.24p we recover the model of the three coupled BEC 
based on the SU{3) symmetry, proposed by Milburn et al. [44] also in the absence of the interatomic 
scattering. 



4.2.3 Two-coupled BEC model with different types of atoms 

We consider the elementary monodromy matrices 

(u + Uai a\a2 aa[ \ / u + Ubi b{ 62 Pb\ \ 

alai u + na2 aal and A[2](7f) = 4^1 u + nb2 pbl 
aai aa2 J \ Pbi Pb2 P'^ / 

T{u) = A^^\u + wi) A^'^\u + W2) (4.25) 

^Notice, however, that it corresponds to a simphficd three-weh model, with no atom-atom interaction terms. 
^More specifically, H = t{0) + mi and we have done the change of variables /?i /3 ; /^a ^ /3 ; /32 — > 7//?, together 
with the following identification Q2 = J ', ^ = ; ^ = {a — /?^) ; ^2 — {a — 7^//?^). 
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Then, the transfer matrix reads 

t{u) ~ u{na + Ub) + UaiUbi + na2nb2 + WiUb + W2na + a\bla2bi + alb\aib2 

+n{a{bi + alb2 + blai + bla2) (4.26) 

where Q = a(3. It leads to the Hamiltonian t(0) with conserved quantities Uj = riaj + Ubj, j = 1,2. 
This corresponds to a model of two wells (^4 and B) with rij atoms of type j, j = 1,2. Here Uaj 
(ribj) denotes the number of atoms of type j, j = 1,2 in the well A (B). Basically, this Hamiltonian 
describes the tunneling of atoms of different types (1 and 2) in the two wells {A and B). 

Notice that this Hamiltonian could also be interpreted as describing two wells (^4 and B) with 
two levels (1 and 2) in each well. Particles can tunnel between the wells and levels. The tunneling 
term Q allow particles to tunnel between wells just in the same level. Similar models (four-mode 
Hamiltonians with tilted potentials) have been proposed recently in [48] . In this context we mention 
that multi-mode models are receiving more attention, specially in connection to the creation of a 
quantum computer from neutral atoms [48]. 

4.2.4 Creation/dissociation of a molecule with two conformations 

We consider the elementary monodromy matrices 

I u + uai A\ai aA\\ / (^a\ (3 AI 

h?-\u) = j a\Ai u + Uai aa\ and A^'^\u) = j p a2 u - na2 ^2 '^2 
\ aAi aai j \ (3 A2 a\A2 u — nA2 

T{u) = A^^\u + wi) A^^\u + W2) (4.27) 

where we have used the sign-transposition and dilatation automorphisms on A[^1(m). Then, the 
transfer matrix reads 

t{u) ~ u{nai-na2)-naina2-Wina2 + W2nai-a'^nA2-/3'^nAi 

+a(3 {A\A2 + AIAi) + a {a\al A2 + A\ 0102) + /9 {a\a\ Ai + A\ 0201) (4.28) 

with conserved quantities 6n = na2 — Uai and ritot = 2 {uai + nA2) + + na2- 

The Hamiltonian describes a molecule A, which exists in two conformations (two different stereo- 
chemical forms) Ai and /I2, and is constituted with two atoms (or submolecules) ai and a2- There are 
transitions between the two conformations of the molecule, and there is recombination/dissociation 
between the atoms and the two aspects of the molecule A. In this context, a and [i are related to the 
probabilities to obtain Ai or A2 starting from the atoms (or submolecules) ai and a2- The relative 
proportion of atoms 5n and the total number of atoms in the system ritot are conserved. 

4.2.5 A generalized heteroatomic-molecular BEC model 

We consider the elementary monodromy matrices 

a2 a 



A'^1(m) = I a\A u + na2 ota\ | ^^"^ ^ 
a A aa2 o? 




T{u) = A^^\u + wi) A^'^\u + W2) (4.29) 
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Then, the transfer matrix reads 

t{u) ~ u{nA + Uai) + ua rial + wi Hal + "-A + ^62 + + P {A^ aia2 + ala{A) 

+ap (4^2 + 4^2) (4.30) 

It leads to the Hamiltonian t(0) with conserved quantity ni = ua + riai- This model describes a 
molecule A which can be dissociated into two different atoms ai and 02- One of these atoms, 02, can 
evolve to a different state 62, which forbids the recombination to A. For instance. A, ai and 02 can 
be trapped in one well, the transitions between 02 and 62 corresponding to a 'leak' of the a2 atom 
toward a second well. 

In the limit — > we recover the model fl4.1.4p . 

4.2.6 Coupling of two oscillators with gl{3) 

For completeness we give here an example of coupling gl{3) with oscillators. We take the case of two 
oscillators, but one can also couple a single oscillator, or three of them. 



Aiu) 



u + Ui a\a2 a a\ 
a\ai 14 + 77-2 a a\ 
aai aa2 o? 



u + en ei2 613 

L(u) = { 621 U + 622 623 

631 632 U + 633 



T{u) = A{u)L{u) 
The transfer matrix reads 



t{u) ~ w (ni + n2 + en + 622) + rii en + n2 622 + a\ 02 621 + 4 «■! ^12 + «^ ^33 
+a{a\ 631 

The conserved quantities are 



+a{a\ 631 + al 632) + a(ai 613 + 02 623) + wi{eii + 622) + W2{ni + 722) 



(4.31) 
(4.32) 

(4.33) 



4.3 Coupling of one oscillator with gl{N) 
We take 

^2 a 



A(7X) 

with 

T(7.) 



/ 



/ ti + en ei2 

621 



eiAT 



(4.34) 



= ftj and Ejj = Cij 



(4.35) 



: A(n)L(7/) 
Then, for /i = 1 and a's real, one gets 

N 

t{u) ~ u (en + ni) + ni en + ^ {cij (a^ e^i + eij a) + a| e^j} + ^ aj ak{ejk + e^j) (4.36) 

j=2 2<j<k<N 

We have not used the parameters Wj, but from the form of A(m) and L{u), it is easy to see that one 
recovers them through the shifts rii — » rii + Wi and ejj ^ Cjj + W2, j = 1,2, . . . , N. 
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4.4 Two-coupled BEC model with {N — 1) levels 

We consider T{u) = AW(m + wi) A^'^^u + W2) with 
^fk = ^jk u + a'^j ttk j,k < N 

= «A^«i ; aSJj-(m) = ftAra] j<N and ASi'^(M) = a 

AfiW = f^^b, ■ A^SA^) = l3^b] j<N and A^S^{u) = (3^^ 
where bj,b^j is another set of oscillators. The transfer matrix t{u) = trT{u) now reads 

N~l 

+ ^ a]akblbj 
where n^^' = aj and n^^' = 6^ 5^. It leads to the 'fundamental' Hamiltonians 

i=i 



m 

t'(0) 



(4.37) 
(4.38) 
(4.39) 
(4.40) 



(4.41) 



(4.42) 
(4.43) 



Again, the shifts nf n'p + Wp, j = 1,2, N - I, p = 1, 2 give back the w dependence. 

It corresponds to a generalization of the previous Hamiltonian fl4.26p , which could be interpreted 
for example, as describing two wells with — 1 levels in each well. The quantities Ij 
j = 1,2, . . . , N — 1, are conserved. 



[1] , [2] 



4.5 iV-coupled BEC model 

We consider the elementary monodromy matrices 



aS(-) 


= 6jkU + ajal 


3,k 


= 1,2 








(4.44) 




= akttj ; A™. 






k 


= 3,.. 


.,N 


(4.45) 


aSH 


= ak aj j, k = 


3,... 










(4.46) 


aS(-) 


= ak aj j, k = 


1,2 










(4.47) 


aS(.) 


= 6jk w + a] ttfc 


j,k 


= 3,...,N 








(4.48) 


aS(-) 


= ajUk ; A|fj 




= "i4 i = li2; 


k 


= 3,.. 


.,N 


(4.49) 



The transfer matrix takes the form 



i(u) 



1=1 i=i 



AT 



aiaj {alttj + ajoj 



rii 



i=l 



(4.50) 
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It corresponds to a multi-well system, which generalizes the previous Hamiltonian (14.241) . The total 
number operator '^f^i Ui is conserved. 



5 Bethe ansatz equations 

We use the results of [49], applied to the special representations we are studying. This can be done 
when the representations are lowest weight. However, when the representations are reducible, the 
Bethe ansatz will give the eigenvalues only on the irreducible parts. It has to be applied for each 
lowest weights. 

For the spin chain part, associated to the elementary monodromy matrix L(m), one has always 
a unique lowest weight v defining the representation (the spin) one is working with. When dealing 
with oscillators, the Fock vacuum |0 > is a natural lowest weight, but different cases can appear, as 
we shall see. 



5.1 N = 2 case 

For N = 2, the elementary monodromy matrices (14. ip and (14.21) become all triangular when applied 
to the Fock vacuum |0 > and/or to the lowest weight vector Vs'- 

Hu)v. ^ (^;^ /^J..,.4z, £(.)|0>^(^ 0)|0> (5.1) 

Aiu) |0 > = f ] |0 > A{u) |0 >= ( ] |0 > (5.2) 



13' J ' ^ ' ' V ^ , 

Note that C-{u) is proportional to the identity matrix, because the representation is reducible. Indeed 
all the vectors |p >= (03)^ |0 >, p G Z+, are lowest weight vectors for C{u): 

«Jf "i-^) |p> (5.3) 

Then, depending on the model we are studying, we will get as pseudo-vacuum for the monodromy 
matrix, either the total Fock vacuum |0,0 >= |0 > CS)|0 >, or different combinations of the type 
\p,q >= \p > ®\q >, \vs,p >= Vs ® \p >, etc... When several pseudo-vacua are at our disposal, we 
will have to repeat the Bethe ansatz method (described below) for each of the pseudo-vacuum in 
order to get a complete set of eigenvalues for the transfer matrix. 

In all cases, the total monodromy matrix is triangular, and the pseudo-vacuum obeys: 



T^j^{u) \vs,Vr > = Xj{u) \vs,Vr > , j = 1,2 with 



Tj^^iu) \vs,p > = Xj{u) \vs,p > , ^^^^ 



\i{u) = {u + wi + s){u + W2 + r) 

\2{u) = {u + Wi - S){u + W2 - r) 

7 = 12 . ( Xi{u) = {u + wi + s){u + W2+p) 
X2{u) = {u + wi - s){u + W2) 



Tf'^{u)\p,q> = X,{u)\p,q> , { i ^ ^\ with 



7 = 12 ( Xi{u) = {u + Wi+p){u + W2 + q) 

' with <^ 

P. 9 e Z+ 1 X2{u) = {u + wi){u + W2) 
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Ai(m) = {u + wi + s){u + W2) 
^2{u) = a [u + wi — s) 



T§^(u)|?;„0> = Xj{u)\vs,0>, j = l,2 with 



r^V)b,o> = A,(.)b,o>, with ; \ ' 

P ^ ^+ X2(u) = a^(u + wi) 



Ai(-u) = (n + u;i)('u + U72) 
' X2{u) = a^f3^ 



T^Mu) |0, > = \j{u) |0, > , j = 1, 2 with 



It imphes that the pseudo-vacuum is an eigenvector of the transfer matrix. We will summarize these 
different cases in the notation 

t{u) \VL >= (^Ai(m) + A2(m)) \n > (5.4) 

where |fi > is one of the pseudo-vacuum(s). 

The other eigenvalues of the transfer matrix then read 

X{u) = Ao{u) Ai(n) + AaH 

where the dressing functions Aj{u) are given by: 

A/ _ 1 _ _ 3 

Aoiu) = Y[l_!^ and Mu) = H^^-^^ (5.5) 

The parameters 1 < n < M, are the Bethe roots, their number M being also a free parameter. 
All these parameters are determined by the Bethe equations 

f^^n-u^-l^X, K + l) i<n<M (5.6) 



5.2 N = 3 case 

The existence of a pseudo-vacuum is not ensured anymore. For instance in the example treated in 
section 14.2.21 one easily computes the action of t{u) on the Fock vacuum: 

t{u) |0, >= {u^ + (3lu + I3l afj |0, > +131 (3; |0, > 

so that the vacuum is an eigenvector of the transfer matrix only if the condition (3i (32C(3 = is 
satisfied. Note that this condition is weaker than demanding |0, > to be a lowest weight vector for 
the monodromy matrix, which would lead to as = 0. 



5.2.1 Pseudo- vacua 

We give the eigenvalues Xj{u) of the pseudo- vacuum (when it exists) under the generators Tjj{u), 
j = 1,2,3. We order them according to the section where they are presented. The pseudo- vacua are 
built on the Fock vacuum 10 >= 10 > ® . . . ® |0 > of the models. 
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Model 14.2. It One must impose the condition (3i(32 = 0. In that case, there are several pseudo- 
vacua \p >= {b^)p |o, 0, >, p e z+. 

Tu{u)\p>=u^\p> ■ T22{u) \p >= Pf{u + p) \p > ■ Ts3{u)\p>=alPl\p> (5.7) 

Model I4.2.2t One must impose the condition /?i/52 = 0. In that case, the Fock vacuum has 
eigenvalues 

Tuiu) \0 >= I3l \0 > ; T22iu) \0 >= piu^ \0 > ; T33(m) |0 >= /^g^ |o > (5.8) 

Model I4.2.3t The Fock vacuum has eigenvalues 

Tn(ii) |0 >= |0 > ; T22{u) \0 >= u"^ \0 > ; T^siu) \0 >= a"^ P"^ \0 > (5.9) 

Model I4.2.4t There are two types of pseudo-vacua 

|p>= (4)P|0, 0,0,0 > and | - p >= (61)^ |0, 0, 0, > , p E Z+ 
Gathering them in the notation |p >, p G Z, their eigenvalues read 
Tii{u) \p >= (3"^ u\p > ; T22{u)\p>=u{u + p)\p> ; Tss{u) \p >= a'^ u\p > , p e Z (5.10) 

Model I4.2.5t One must impose the condition aP = 0. In that case, the Fock vacuum has eigen- 
values 

Tii(m) |0 >= |0 > ; T22{u)\0>= P'^u\0> ; Tss^u) \0 >= a'^ u\0 > (5.11) 

Model I4.2.6t We denote by v the 5^/(3) lowest weight vector, with eigenvalues (Ai, A2, A3) under 
the gl{3) Cartan generators. The pseudo-vacuum is \v >= \0 > ^v, with eigenvalues 

Tn{u) \v >= u{u + Xi) \v > ; T22{u) \v >= u{u + X2) \v > (5.12) 
n^iu) \v >= a^u + X3) \v > (5.13) 

5.2.2 Bethe equations 

Now that we have determined the conditions for the existence of a pseudo-vacuum, one can write 

t{u) \n >= (^Ai(w) + A2(w) + A3(ii)) \n > (5.14) 

where the eigenvalues Xj {u) can be read in the section 15.2.11 Then, the other transfer matrix eigen- 
values read 

X{u) = Ao{u) Xi{u) + Ai{u) X2{u) + A2{u) X-siu) 
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with the dressing functions 



mW (1) , 1 



n=l ^ u,n 2 

A/(i) (1) 3 Af(2) 
M-M^ -^ T-r U-Ur' 



U Un 2 n=l U Un i 

^''^ _ ,,(2) _ ^ 

= n I) • (5-17) 

„=1 M - ^ - 1 

We have here two sets of Bethe roots, u^n\l<n< M(^) and Un \ M^'^\ constrained by the 

Bethe equations 

TT TT = — ) 4- l<n< M« (5.1?^ 

J-A ^,(1) _ ,,(1) I 1 J-i (i)_ (2)_i > / ai . A - - ^ 

^(2) _ (1) 1 (2) (2) _ 

n ; ; n ;;, = , , i<n<M<^) (s^w) 

5.3 General case 

For the general case of N x N matrices, and supposing the existence of pseudo- vacua 1^2 > (possibly 
with conditions on the parameters of the models, as above), we will get for the eigenvalues of the 
transfer matrix 

N 




k=l 



with pseudo-vacuum eigenvalues 



Tkkiu) \n >= Xk{u) \n> , k = l,...,N so that t{u) \n >= Xk{u)^ \^ > 



and dressing functions 



Mn)=u I n o<k<N-i, (5.20 



n=l U Un 2 n=l U Un ^ 



with = M(^) = 0. 

The — 1 types of Bethe roots u^n^ , 1 < n < M^^\ \ <k < N - \ will be determined by the 
Bethe equations: 

MC^-D (fc)_ (fc-l) 1 Ak) W _ J'^+l) , 1 \k{u^n 

7w TTfc^^^i) T 11 7(fc) „ (fc) , 1 11 ~ 



J-l (fc)_ (fc) . 11 , r^W, fc^ 

1 < n < M^*^) and 1 < < iV - 1 (5.21) 
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They depend on the model and pseudo- vacuum through the eigenvahies Xk{u), k = 1, . . . , N. 

6 Superalgebras and fermions 

We have seen that the Yangian is not sufficient to build BEC models based on fermions. For such 
a purpose, one needs to consider another algebraic structure, the super- Yangian Y{M\N) based on 
the superalgebra gl{M\N). 

The models look very similar to the ones presented in section HI with the notable difference that 
the choices for the fermionic oscillators is not the same, allowing more flexibility in the construction. 

6.1 The superalgebra gl{M\N) 

To define a superalgebra, one needs to introduce a grading [.] that will distinguish the fermionic 
generators from the bosonic ones. For gl{M\N), denoting Cij, i,j = 1, ...,M + N these generators, 
the grading is defined a^ 

,,,,,, , fo if 1<A;<M , , 

k, = [i] + hi with [k] = { - - 6.1 

The generators with grading are bosonic ones; they form a gl{M) © gl{N) subalgebra of gl{M\N), 
generated by Cij, with i,j < M and i,j> M. The remaining generators are of fermionic type, and 
form a representation (A^, M) © (A^, M) of this subalgebra. The supercommutator reads 

[e., , e,,} = -(-l)(W+b1)(W+m) [e,, , e,,} = 5,, e, - (-l)(W+b1)(W+H) 

which amounts to consider commutators when e^j and/or Cki are/is bosonic (i.e. of grade 0), and 
anti-commutators when both are fermionic. 

The fundamental representation is of dimension M + N: 

Ti{eij) = Eij 

where now Eij are graded elementary matrices of size M + N. 

6.2 The super- Yangian Y{M\N) 

One defines the super- Yangian through a graded /2-matrix, obeying a graded version of YBE. By 
graded version of YBE, we mean that one has to use a graded tensor product on the auxiliary spaces: 

{E,, © Eki) ■ {Ep, © Ers) = (-1)(W+W)(W+M) {E,,Ep,) © (EkiErs) (6.2) 

where the grading is the same as the one of gl{M\N). 

In fact, everything looks formally the same as for the Yangian, with the restriction that one has 
to take care of the grading (16.11) . For instance, the i?-matrix reads: 

i?i2(x) = I©I- -Pi2, (6.3) 
X 

"other choices of grading could be used, but we wiU stick to this one throughout the paper. 
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with now the super-permutation operator 

M+N M+N / M M+N 

i,j=l i=l \j=l j=M+l 

Plugging this i?-matrix in the relation (12.11) . and taking care of the graded tensor product leads to 

[Li,{u) , Lki{v)} = ^ ' [L^,j{u)Lu{v) - Lkj{v)Lu{u)j . (6.4) 

We will use elementary monodromy matrices built on gl{M\N): 

Lij{u) = u5ij + Cij , Cij e gl{M\N) (6.5) 

6.3 Transfer matrix and symmetries 

The monodromy matrix /S.^^'>T{u) = T^^^ (u) . . .T^^^ (u) gives a transfer matrix of the form 

st{u) = str{T^^\u)...T^^\u)) (6.6) 
where the super-trace of a matrix is defined by 

M+N M M+N M+N 

str(A)= J](-1)[^U,, = ^A,,- ^=E^^'^-^*^- (6.7) 

1=1 i=l i=M+l «J=1 

As for the bosonic case, one can show that 

[st{u) , st{v)] = (6.8) 

Again, to get Hermitian Hamiltonian, we will focus on the case L = 2, so that the monodromy and 
transfer matrices we will be concerned of, have the form 

2 

Tm{v) = ^t;"r,f) with Tl^f=uJkSki and ook e C (6.9) 

n=0 

st{u) = t2u'^ + tiu + tQ with t2 e C (6.10) 

Specializing to gl{M\N) representations will give different models. These models will have also a 
symmetry, as one can see from the relation 

[st{u) , ri;)] = {u, - ui) (Ti?) + «ri;)) (6.11) 

proving again that the quantities 

h = Tii , VA; = 1,...,A^ (6.12) 

commute with the transfer matrix. 
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We will be essentially interested in the oscillator representation 

7i{eij) = a] ttj with [a^ , a]} = /ij^jj (6.13) 

where (aj,a|) for 1 < i < M (resp. AI+l<i<M + N) are bosons (resp. fermions), i.e. 
h] = [al] = [^]. 

This choice of grading implies that [eij] = [a]] + [aj] = [i] + [j], in accordance with the gradation 
of gl{M\N). It corresponds to a bosonic (resp. fermionic) oscillator representation for the gl{M) 
(resp. gl{N)) subalgebra. Remark that the opposite choice (i.e. [aj] = [aj] = [i] + 1) is also possible. 



The inhomogeneous oscillator monodromy matrices will then be obtained by taking constant 
bosonic oscillators. In that process, the choice of the grading for the oscillators will be essential, 
since it will determine which of the oscillators can be possibly set to constant. 

Despite of the grading, the elementary matrices are still hermitian. Hence, the models will be 
hermitian. 



7 Fermionic BEC models 



7.1 Two by two (graded) matrices 



We are dealing with the gl{l\l) case, and we focus on the elementary monodromy matrices of the 
form 



JZiu) 
A{u) 



U + TLi 


a\ 02 




u-n2 



u + n 


/3ct ~ 


(3c 





A{u) 





Pc 




u — n 



(7.1) 
(7.2) 



where in C{u) one of the couple (aj,aj), j = 1,2, is bosonic and the other one fermionic, while in 
A(u) and A(w), (c, c^) is fermionic. 



7.1.1 Fermionic heteroatomic BEC model 

We consider strC^^Mu + wi) 

£[2] 

{u + W2), with oscillators (ai,ci) for £W(m) and (02,02) for C^^\u). 
One can consider either (ai, 02) to be bosonic and (ci, C2) fermionic, or (oi, 02) fermionic and (ci, C2) 
bosonic. 

The transfer matrix reads: 

St{u) ~ U{nal + na2 + rtcl + nc2) + r^al rtcl - na2 nc2 + Wi(ncl + 71^2) + «'2('^al + na2) 

+a\clcia2 + alc\c2ai (7.3) 

We get a fermionic version of the model described in section 11.1.31 It can be interpreted as modelizing 
a coupled pair of one boson and one fermion which can tunnel together from one well to another. 
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7.1.2 Heteroatomic-molecular BEC model with fermions 

We consider strC{u + Wi) K{u + W2), with oscillators (ai = b,a2 = a) for C{u) and c for A(m). In 
jC{u), we take 6 bosonic. Then, a is fermionic, as well as is c in A{u). We get 

st(ti) ~ u{nb + ric) + nbUc + wiHc + W2nb + fS"^ Ha + f3 {a^c^b — b^ac) (7.4) 

We find a new version of the heteroatomic-molecular BEC model of section 14.1. 41 with now a bosonic 
molecule b constituted of two fermionic atoms a and c. 

The shift automorphism applied on b produces additional terms in the transfer matrix 

Hhound ~ {P{b + b^) + (3^)n, + W2f3{b + b^) + a(3{ah^ -ac) (7.5) 

7.1.3 Fermionic two- wells 

We start with A^^^(u + wi) and A^'^^{u + W2), set all the bosons to constant, keeping the fermions ci 
and C2- 

strA^^\u + wi)A^^\u + W2) = iu + wi + nci){u + W2 + 12^2) - Pip2c\c2 + P1P2C1CI (7.6) 



st{u) ~ u{nci + nc2) + ncinc2 + Winc2 + W2nci- PiP2{clc2 + clci) (7.7) 

One recognizes a two-wells models, but now with atoms of fermionic nature. 

We mention here that there has been recently a great interest in Bose-Einstein condensates with 
fermions since the achievement of quantum degeneracy in ultracold Fermi gases (see, for example [18] 
and references therein). 



7.2 Three by three (graded) matrices 

We are now dealing with the 5f/(2|l) case, corresponding to 3 x 3 matrices. Again the number of 
possible models becomes numerous, so that we present the generic case and treat only two examples, 
physically relevant. 



7.2.1 Generic case 



We present here a general formulation for the transfer matrix, which encompasses all the possible 
models, by setting some of the (bosonic) oscillators to constant. 
The two generic elementary monodromy matrices have the form 



lW(«)= 1 




h ni 
ai 


a\ a2 
U2 + n2 


a[ 
al 


as ^ 

as 


and LP](n) = | 


I' vi+rrii c\ C2 
c\ Cl V2 + m2 


c\cs 

4 C3 




K 4 




a\ a2 


Us H 


- ns / 




\ C3 Cl C3 C2 


Vs + mg 


with the notations 






a] aj 


and 


rrij = c] Cj j = 


1,2,3. 





(7.8) 

For LW(m) to be of ^/(2|1) type, one can choose either ai and 02 bosonic and 03 fermionic, or a\ and 
a2 fermionic and as bosonic. Obviously, the same criteria apply for LP\u), Cj, j = 1,2,3. 
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When the oscillators are bosonic, one can choose to set them to constant: in that case the 
corresponding spectral parameter {uj or vj) has to be set to zero. In all other cases, it is set to 
u. For instance, if oi, bosonic, is the only one set to a constant ai, then, one has ui = and 
U2 = = vi = V2 = = u. We will also have in this case a\ = al and rii = 

Keeping these rules in mind, one can compute a generic transfer matrix. It takes the form (up 
to irrelevant constant terms) 

3 

st{u) = ''^^{uj nij + Vj rij) +H 

3 

H = a\a2 c\ c\ + c| C2 a\ + a| as Cg c\ + c| C3 Og a\ + as Cg C2 + 4 C3 Og a2 + rij 

j=i 

The 'true' Hamiltonian st(0) = if of a given model is then obtained through the above rules, after 
choosing which of the oscillators are bosonic or fermionic, and, among the bosonic ones, which of 
them are set to constant. 

The conserved quantities of the model take the form 

u Ij = Uj rrij + Vj rij , j = 1, 2, 3 

Of course, depending of the choices, some of these quantities can be zero after use of the rules. 

Let us also remark that if one takes all the oscillators to be bosonic or fermionic (a situation 
forbidden in the case of gl{2\l)), one gets a generic transfer matrix for the usual three by three 
matrices, as treated in section 14. 2[ 

7.2.2 A fermionic two-coupled BEC 

To exemplify the above techniques, we construct a fermionic version of the model given in section 
14:2:31 

We take as and cs to be bosonic, and set both of them to constant, a and 7 (both real) respectively. 
Then, one gets four couples of oscillators, all of them being fermionic. The rules lead to 

Ui = U2 = Vi = V2 = u and = v-^ = Q (7.9) 

as = ag = a G M and = o? (7-10) 

Cs = 4 = 7 e M and ms = -/"^ (7.11) 

Thus, we get an Hamiltonian 

H = a[a2 c\ ci + c\ C2 al ai + 07 {a[ ci + c\ai + ^2 + 4 a2) + mi ni + 1112 n2 
with conserved quantities 

11=1^1+ mi and /2 = ri2 + m2 

We recover the Hamiltonian of section 14.2.31 with the notable difference that the atoms have a 
fermionic nature. 
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7.2.3 Fermionic three coupled BEC model 

If now one takes 03 bosonic and constant (a G M), and C\ and C2 bosonic and constant (7j G M), one 
gets: 

ux= U2 = v-i = u and u->, = = V2 = ^ (7-12) 
as = 03 = a G M and 713 = (7-13) 
Cj = c] = G M and rrij = 7^ , j = 1, 2 (7.14) 

leading to (using the sign-transposition and dilatation automorphisms on //'^'(w)) 

H = 7i72 (al a2 + 4 «i) + "7i («i C3 + 4 «i) + "72 (4 C3 + 4 02) + 7i ni + 72 n2 - a ms (7.15) 
We get the fermionic version of the Hamiltonian described in section I4.2.2[ 

8 Conclusion 

In the BEC context, we have constructed integrable generalised models in a systematic way exploring 
different representations of the gl{N) algebra and the gl{M\N) superalgebra. Some existing models, 
such as the two-site Bose-Hubbard model, have been recovered and new ones have been predicted. 
Interestingly, a two-coupled BEC model with a field, a three-coupled BEC model and a two-coupled 
BEC-model with different types of atoms, among others, have been introduced. The use of the 
gl{M\N) superalgebra allows the introduction of fermions, leading to systems mixing bosons and 
fermions, as they are presently studied in condensed matter BEC experiments. In this context, the 
'integrable approach' can be viewed as a technics to construct in a very general way Hamiltonians 
relevant for these studies. 

The energy spectrum of these models has been derived, through the Bethe ansatz equations, 
by the use of analytical Bethe ansatz. The next step in the study of these systems in this general 
framework, is the determination of the (Bethe) eigenstates and eigenfunctions, which would allow to 
investigate the classical and quantum dynamics of such systems. 

Finally, we remark that more general integrable models can be obtained using the method pre- 
sented in the present work. They are constructed using products of more elementary monodromy 
matrices, with the restriction that the hermiticity of their Hamiltonian is not guaranteed anymore. 
Apart from the trial and error method that one can use on a case-by-case basis, a general analysis 
determining the conditions under which Hamiltonians are hermitian would certainly improve the 
landscape of integrable BEC models. 
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A Quantum determinant and conserved quantities 



We have seen that the expansion of the monodromy matrix provides some conserved quantities of 
the integrable models. However, other conserved quantities can be obtained when considering the 
center of the algebra. For the Yangian, it is known that its center is generated by the quantum 
determinant [50]: 

qdet(u} = ^ sgn{a) Ti,^(i)(ti) T2^a{2){u - 1) ■ ■■TN,a{N)iu - + 1) (A.l) 

where Sn is the group of permutations. It is clear that the conserved quantities obtained in this 
way are quite complicated, but they may be of some help for the study of the different models. To 
illustrate this, we give the form of the quantum determinant when N = 2 and 3. 

qdetiu) = Tn{u)T22{u-l) -Tu{u)T2i{u-l) for = 2 (A.2) 
qdet{u) = Tn(ti) T22(n - 1) T^siu - 2) + T^2{u) T23{u - 1) Tsi{u - 2) + Tri{u) T2i{u - 1) T^2{u - 2) 
- T12W) T2i{u - 1) T^^{u - 2) - Ti3(m) T22{u - 1) T32(^i - 2) - rn(M) T2'i{u - 1) T;2{u - 2) 
for A^ = 3 (A.3) 

qdet{u) is a polynomial in m, of degree 2A^ since T{u) is of degree two. Hence, one gets a priori 2A^ 
conserved quantities. They are not all independent, but they can provide new conserved quantities, 
not contained in the transfer matrix, nor given by (13.61) . Remark that these quantities are a priori 
not hermitian, but, since they are central in the whole Yangian, so are their adjoint. Hence, one can 
build hermitian (and anti-hermitian) conserved quantities from qdet{u). 

Indeed, for N = 2, expanding qdet{u) from the expansion (13.21) . one gets (up to constant terms) 

3 

qdet{u) = (A.4) 

n=0 

ds = u^Ti^+U2Til^ (A.5) 

d2 = u.T^^ -2uj2Til^ +u;^Ti^ +U2Tif -T^^tS^ (A.6) 

d, = U2 - 2.2 Ti? + Til^ ri? + T,?) - Tg) T^? - T^^ (A.7) 

do = a;2r«-TS)r« + T(?T« (A.8) 

with Uk = fiak fJ'bk, k = 1,2, wheu dealing with oscillator representations, or cj^ = 1 for gl{N) 
representations. 

After some algebras, and using the conserved quantities (13. 6p . one obtains the following invariants 



Co 


= Ti^^Til^ + TS\u2 


-r«) 


(A.9) 


Cl 


= T,?T2?) + r(?(-i 


-T«) 


(A.IO) 


C2 


^(1) ^(1) ^(0) 


-U.2TS' 


(A.ll) 



Of course, the explicit form of these invariant will depend on the representations we will use, i.e. on 
the physical model we are studying. 



28 



In the case = 3, the same kind of calculation leads to more complicated expressions. The 
simplest ones read: 

C4 — ^3 -'12 -'21 +^2-^13 -'Sl +'^l-'23 -'32 —^2 1^3-' 11 — ^1 l^3 22 — 1^1 Ci^2 ^ 33 l^-J-^J 
^ _ ^(0) ^(0) ^(0) _ ^(0) ^(0) ^(0) _ ^(0) ^(0) ^(0) ^(0) ^(0) ^(0) ^(0) ^(0) ^(0) _ ^(0) ^(0) ^(0) 
t-O — 11 -'22 33 11 -'23 -'32 12 21 33 "T 12 23 31 "T 13 21 32 -'13 -'22 31 
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